Phase separation in effective hard-core boson and triplet models in one and two 

dimensions 
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Effective models of hard-core hole pair bosons and triplets are derived from the t-J model in ladders 
and on the square lattice by performing a change from site to dimer basis. The Hilbert space is 
truncated by projecting out single-occupied electron states and only nearest neighbor interactions 
are retained. The resulting effective models in one and two dimensions are studied by numerical 
techniques. In both spatial dimensions, the main result is that each hole pair is surrounded by a 
singlet cloud expelling triplet excitations from its vicinity. It is suggested an interpretation of this 
feature as a phase separated state between a hole-pair rich singlet phase and an undoped triplet 
phase with antiferromagnetic correlations. The possible relevance of this result to other theoretical 
scenarios and experimental results is discussed. 

PACS: 74.20.-z, 74.25.Dw, 74.25.Ha, 74.20.Mn 



I. INTRODUCTION 

The S0(5) theory! is a very appealing framework in 
which antiferromagnetic (AF) and superconducting (SC) 
phases of high-Tc superconductors are naturally related. 
This theory not only gives an exulanation of the resonant 
peak observed in the SC phaseo but also predicts that 
it smoothly connects with the magnon peak of the AF 
phase as it was recently observed experimentally.u An- 
other of the predictions of this theory, i.e. the presence 
of AF order inside the vortices present when a magnetic 
field is applied on the SC phase_has recently motivated a 
flurry of experimental activityu'El confirming partially this 
prediction. At the realm of the SO (5) theory lies the fact 
of singlet pairing of electrons. This is a deeply strongly 
correlated electrons feature which shares with resonant 
valence bond (RVB) theoriea3 and excludes an explana- 
tion of the most important features of SC cuprates, not 
only the AF and SC phases but also the intervening pseu- 
dogap phase, via Fermi liquid or Fermi liquid instability 
concepts. Although there have been many attemp^of 
formulating a S0(5) symmetric model on lattices,Lra a 
connection between the continuum theory and a more 
microscopic model like the t — J model is still missing. 

In this sense, one of the motivations of the present 
study is to help bridge the gap between the more phe- 
nomenological S0(5) theory and the microscopic t — J 
model. More specifically, our goal is to obtain and study 
an effective model obtained by a change of basis from 
the site, spin-1/2 electrons, basis to the dimer basis, and 
then projecting out single occupied dimer states. This 
dimer basis also indicates a relation with RVB. scenarios, 
in particular with its nearest neighbor version.Q However, 
in the present work the change of basis is performed on 
a single spatial dimer covering and hence the problem 
of an overcomplete basis which affects RVB models is 
avoided. In addition, the elimination of single-occupied 
dimers implies that the excitations of our model will not 



be spinons and holons, resulting from the fractionaliza- 
tion of electrons, but, as in the S0(5) theory, collective 
states of singlet binding of electrons in dimers id As a 
matter of fact, a recent study of the spin-1/2 Heisenberg 
model in two dimensions (2D), starting from a singlet 
RVB "soup" , has precisely arrived at the conclusion that 
the elementary excitations of the AF gmund state are 
spin-1 excited dimers instead of spinons Jlj 

The pseudogap phaseO, to many a central piece of the 
puzzle of understanding high-Tc superconductivity, is a 
candidate to be a realization of a SO (5) symmetric state 
or an RVB state. An effective model on dimers could then 
be a good starting point to understand many featureSpOf 
the pseudogap phase. In addition, the stripe phaseEJ, 
which in some cuprates appears inside the pseudogap 
phase, could be explained in terms of a particular or- 
dering of dimers that is pinned by a spa. tiall jy anisotropic 
exchange or by coupling to the latticcll^llil In fact, the 
concept of electron singlets plav-a central role in one of 
the theories of the stripe phase.Ej 

There is a more specific, recent experimental result 
which could be possibly described in terms of an effec- 
tive model oir dimers. A scanning tunneling microscopy 
study of Bi2Sr2CaCu208+5 has shown a phase separa- 
tion between superconducting islands inside a percolat- 
ing background with phjsical properties resembling those 
of the pseudogap phase. 

Other motivations to study effective models are both 
technical and conceptual. From the methodological point 
of view the reduction of the Hilbert space is an advan- 
tage to most numerical techniques, exact diagonaliza- 
tion, density-matrix renormalization-group, and quan- 
tum Monte Carlo. Conceptually, the effective model may, 
as it is shown below, make evident some properties which 
are somewhat hidden in the original t-J Hamiltonian. 

As said above, our approach to SO (5) concepts is to 
map the 2D t-J model onto an effective model in a basis of 
singlet electron dimers. Since the original model excludes 
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double-occupied sites, this connection is actually closer to 
a recently proposed variant of the S0(5|)-i»©del called the 
"projected" S0(5) (or pS0(5)) modeLO'ta This pS0(5) 
model is an attempt to close the bridge with a micro- 
scopic model on a lattice from the other end with respect 
to what is intended in the present work. Although the 
basis set of the pS0(5) model is identical to that of the 
models that will be shown below, there are important 
differences in the corresponding Hamiltonians. 

The paper is organized as follows. In Section ||, an 
effective model is derived for the t-J model on ladders. 
In this case, the procedure can be carried on in a much 
cleaner way than for the t-J model on the square lattice. 
In addition, the one-dimensional (ID) character of this 
effective model allows its study by exact diagonalization 
techniques. The main results are common to the corre- 
sponding ones for the most interesting case, i.e. on the 
square lattice, which are obtained in Section III . Finally, 



in the Conclusions, the possible relevance of the present 
results to understand a number of previous theoretical 
problems as well as experimental results is discussed. 



II. LADDERS 

Although the main interest in connection with high- 
Tc superconductivity is to study the interplay between 
various strongly correlated electron phases in 2D, it iMn- 
structive to start with the analysis of two-leg ladders.Ej It 
should be noticed that the t — J model on this lattice has 
been initially studied, among various reasons, as a real- 
ization of an RVB state.EB In addition, extensive studies 
have revealed strong similarities between the behaviors 
obtained on ladders and on the 2D square lattice. 

The Hamiltonian of the t — J model is: 



+ Jl_Y,i^t,^~l) (2) 

i 

where tj^, ^ is a creation operator of a triplet with — 

0, 1, —1 {a = 0, +, — ) at dimer i, nn — riQj + ri_|-.i + 

^ci = t]^ ita,i- The first term corresponds to the spon- 
taneous creation of two triplets out of the vacuum with 
the constraint of keeping total S*^ = and the second 
term is a triplet "hopping" . It has been noted that this, 
constraint is not satisfied in the pS0(5) Hamiltonian.E3 
The third term is the usual spin- 1 Heisenberg exchange 
Hamiltonian. This Hamiltonian (or equivalent expres- 
sions) has been derived previously (e.g. Ref. ||, and 
references therein), and used in restricted Hilbert space 
diagonalizations.Eil The Hilbert space of the effective 
model turns out to be disconnected into two pieces: the 
set with even and the set of odd number of triplets. For 
example, in the subspace of total — the first set 
is generated by successive applications of the Hamilto- 
nian to the initial state {. . . sssss . . and the second 
one to the initial state {. . . sst^ss . . where s stands 
for a singlet dimer. This splitting of the Hilbert space 
does not appear in the effective model obtained for chains 
(Fig. |l|(a)). 

Away from half-filling, the states with a single hole 
occupancy on a dimer are projected out. If V is the 
projection operator on the subspace of retained states 
and Q is the projection operator on the subspace of the 
eliminated states, then the effective Hamiltonian is given 
by the standard formula: 
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n^Tij) (1) 



where the notation is standard. On ladders, = t, Jij = 
J along the legs, and tij = Jij = J± on the rungs. 
On chains and on the square lattice, the isotropic and 
homogeneous case {tij — t, Jij — J) will be considered. 
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FIG. 1. Change from a site to a dimer basis in (a) a chain 
and (b) a ladder. 

At half-filling, the exact change of basis from the site to 
the dimer basis (Fig. |l|(b)) leads to the following Hamil- 
tonian: 
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J.ladder 
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where H^-o = £^o*o, and HeffV-^o = ^^o'PI'o- If H is 
the one given by Eq. (0) on the 2 x 2 ladder with 2 holes, 
the effective Hamiltonian will contain nearest neighbor 
(NN) terms only. In this way, the hopping term of the 
effective Hamiltonian, 'Ht,eff results: 



t^ladder 



tsY^iblh + h.c. + 

J±J2(np,,-l) (4) 



where fej is a creation operator of a hole pair in the dimer 



= b]b. 



— 1,0 if the site is empty (i.e. 



singlet) or occupied (by a triplet or a pair) respectively. 
The first term corresponds to the hopping of a pair to 
a singlet site, while the second term corresponds to the 
hopping to a site occupied by a triplet. The coupling 
constants ts,tt are complicated functions of the original 
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parameters {J, J±,t,t±} but always tg > tt, as shown in 
Fig. ^ for two values of the lattice anisotropy a = t±/t 
{ts — 2t(, for the isotropic ladder a — 1). This is true 
even in the presence of a NN pair-pair Coulomb repulsion. 
This relation between the hopping parameters already 
suggests one of the most important results of this work. 
That is, the pairs would tend to be surrounded by singlets 
rather than by triplets in order to gain kinetic energy. In 
addition, the Heisenberg term in Eq. would favor the 
clustering of triplets. These two combined effects would 
imply the phase separation between a pair-doped singlet 
phase and an undoped antiferromagnetics phase. 

Finally, the full effective Hamiltonian in the projected 
dimer basis is given by: 

-p/ieff) _ njidimer) ^{eff) ,r\ 
''■ladder — ''■J,ladder '^tdadder 

The ID effective model can be studied by exact di- 
agonalization techniques (Lanczos algorithm). Most of 
the results below were obtained on a L = 12 chain with 
periodic boundary conditions. Some computations for a 
L = 16 chain show that finite size effects are not impor- 
tant. All results shown below correspond to fixed num- 
ber of pairs Np, and were obtained for the ground state, 
k = (0,0). In the following, all energies and coupling 
constants are expressed in units of t of the original t-J 
Hamiltonian Eq. (|l|). 

2.4 I ' • ' • ' • ' • ' 1 

sp <-> ps (a=1) 

2.0 tp<->pt(a=1) 

sp <-> ps {a=2) 

tp <-> pt (a=2) 
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FIG. 2. Hopping amplitudes between pairs p and singlets 
a (island between hole pairs and triplets t {tt), as defined in 
Eq. M), as a function of J and for a = 1 and 2. Basis change 
on a ladder (Fig. |l](b)). 

In Fig. ^ the correlation between a pair and a triplet, 
(ntUp) at distance r — 1 (NN sites) is shown for the 
L = 12 chain with one pair, as a function of J/t and 
for three ladder anisotropy ratios. It can be seen that 
the probability of finding a triplet next to a pair is much 
smaller than the probability of finding a pair next to 
a singlet. The later, in the one pair system, is simply 



(nsUp) = 1 — (ntUp). It is more important the result 
that, as it can be easily seen in Fig. [ntUp) (r = 1) 
is always smaller than the triplet density {nt) (inset) for 
the same J/t. The normalization of the later is such that 
ns+nt + rip — 1, Up — 0.0833 (one pair), and Up — 0.1667 
(two pairs). This result suggests that triplets are expelled 
from the vicinity of a pair. 
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FIG. 3. Triplet-pair correlation at r = 1, L = 12, A'p = 1 
(open symbols), as a function of J and for the values of the 
anisotropy ratio a indicated on the plot. The inset shows the 
triplet density for the same cluster and parameters. Results 
for Np = 2 (filled circles), a = l axe also included. 
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FIG. 4. Diagonal correlations as a function of distance ob- 
tained for L = 12, A'p = 1, a = 1, and various values of 
J, as indicated on the plot, (a) Triplet-triplet (dot lines) 
and triplet-pair (dashed lines) correlations; (b) singlet-singlet 
(full symbols, divided by 2) and singlet-triplet (dashed lines) 
correlations. 
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The complete picture can be inferred from Fig. ^ 
by looking at various correlations {uarib), (a, 5 = 
to,t+,t_,p) as a function of distance. From now on, the 
study will be limited to the isotropic case {t ^ t±, J = J± 
in the original t-J model) but similar results were also 
found for all a > 1 investigated. In this Figure, which 
corresponds also to the system with one pair, it can be 
observed that the triplet-triplet correlation is maximum 
at 7' = 0, while the triplet-pair correlation is maximum 
at the largest distance on the chain. That is, triplets try 
to stay as far apart as possible form a pair. On the other 
hand, although singlet-singlet correlations are also maxi- 
mum at r = 0, the singlet-pair correlations are maximum 
at r = 0. The final piece is that the singlet-triplet cor- 
relation is maximum at the maximum distance. Similar 
results were obtained for the larger L = 16 site chain. 

The picture emerging from these correlations is that 
pairs move preferentially in a background of singlets, and 
that both pairs and singlets try to keep themselves away 
from triplets. Thus, the system is separated between a 
pair-rich singlet phase (it is tempting to consider this 
phase as a doped "RVB" phase, even though that the 
present effective models are obtained for a single dimer 
covering) and a pair-poor triplet-rich phase which can be 
identified as an undoped AF phase. 
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FIG. 5. Diagonal correlations as a function of distance ob- 
tained for L = 12, Np = 2, a — 1, and various values of J 
as indicated on the plot, (a) Triplet-triplet (full symbols), 
pair-pair (open symbols) and triplet-pair (dashed lines) cor- 
relations; (b) singlet-singlet (full symbols, divided by 2) and 
singlet-triplet (dashed lines) correlations. 

Essentially the same behavior is observed in the case of 
two pairs present in the system. Results for the triplet- 
pair correlation in NN sites and the triplet density as a 
function of J/t are also included in Fig. ^ for comparison 
with the one pair case. Both the probability of finding a 
pair next to a triplet and the triplet density are smaller 
than for the one pair case. A likely explanation of this 



behavior is that the introduction of more pairs in the 
system increases the volume of the pair-doped "RVB" 
phase leaving less room for the AF phase. 

In Fig. H, diagonal correlations are shown for L = 12 
and same parameters as in Fig. ^. The new feature in 
Fig. ^ with respect to Fig. ^ is the fact that two pairs re- 
pel themselves as can be read from the fact that pair-pair 
correlations are maximum at the largest distance. Then, 
the same qualitative behavior found in the one-pair case 
holds: the system separates into a pair-doped singlet re- 
gion, here formed by two islands and a pair-poor triplet- 
rich region, presumably with short-range AF order, in 
this case filling the space between those two islands. No- 
tice also in Fig. ||(b) the singlet-singlet correlations falling 
down to its bulk value within a lattice spacing. 

III. TWO DIMENSIONS 

The most important situation is that of the square lat- 
tice, which corresponds to the Cu02 planes in super- 
conducting cuprates. There are again infinitely many 
different ways in which a change from the site to the 
dimer basis can be performed. One of them is shown 
in Fig. A well-known feature of choosing a dimer ba- 
sis like the one depicted in Fig. ^ is that the rotation 
invariance of the square lattice is broken and it is not 
simple to restore it at the level of the effective Hamil- 
tonian. In the ladder case, examined in the previous 
Section, this is not important since the lattice itself is 
spatially anisotropic. In the case of the square lattice, 
the purpose of the present study is to provide indications 
of the presence, in a rotationally-broken effective model 
for the square lattice, of the singlet-AF phase separa- 
tion already observed in ladders and to suggest that this 
feature should be also present in a rotational invariant 
formulation£3 




FIG. 6. A possible change from a site to a dimer basis in the 
square lattice. Dotted and dash-dotted lines indicate effective 
magnetic interactions in the new basis. 

At half-filling the change of basis implies an exact map- 
ping of the original spin- 1/2 Heisenberg model onto a 
model written in terms of triplets and singlets. In the 
basis change indicated in Fig. ^ the resulting model is 
anisotropic with ladder-like interactions (Fig. ^b)) given 
by Eq. (H) in the horizontal direction and chain-like in- 
teractions (Fig. |l|(a)) in the vertical direction, which are 
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given by the following Hamiltonian: 



{eff} 



n 



{dimer) 
J,chain 



J±E("M-1) (6) 



where /3 = +, 0, — , 7 = — , +, for a = 0, +, — respec- 
tively. The sign of the second term depends on the defi- 
nition of the singlet. It should be noticed that this term, 
as in the ladder case, locally conserves the total 5^. This 
Hamiltonian can be read, with a slightly different nota- 
tion though, in Ref. |o|, and alternative or similar deriva- 
tions can be found in several other places .ESEil Finally, 
the exchange part of the Hamiltonian in the dimer basis 
is: 



n 



(di', 



r) 



J^square 



n/(dimer) . q_,{dimer) 
Jjladder J,chain 



(7) 



with the ladder (chain) term acting on the horizontal 
(vertical) direction as indicated in Fig. O. 
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FIG. 7. Hopping amplitudes between pairs p and singlets s 
{ts,tso) and between pairs and triplets t {tt,tto), as defined in 
Eq. (|^ , as a function of J. Basis change on a chain (Fig. 0(a) ) . 

Let us consider the original t-J Hamiltonian defined 
on a chain. To obtain the hopping part of the effective 
model, the projection defined by Eq. (^ is again used. 
In this case, H is the t-J Hamiltonian on a 4-site chain 
(Fig. |l|(a)) with two holes and the effective one will con- 
tain NN interactions only. The resulting effective hop- 
ping term is very similar to the one in Eq. (^: 



t^chain 



< E i%A,sh + h.c.) 



+J±^{np,i - 1) 



(8) 



As in the case of Eq. (Q), the hopping amplitudes sat- 
isfy t'g > t[, and tsQ > ttQ (Fig. 0) again favoring the 
movement of pairs away from triplet-rich regions. 

Coming back to the square lattice, the correct pro- 
cedure is to to take Ti as the t-J Hamiltonian on the 
eight-site cluster indicated by a dashed box in Fig. |[ As 
a result, the effective Hamiltonian contains three- and 
four-site terms in addition to NN interactions. In order 
to keep the HamiltoniaiJr-as, close as possible with the 
proposed pS0(5) model,ll3tJ only NN hopping interac- 
tions as given by Eqs. (^) and (^) in the horizontal and 
vertical directions respectively are retained. The signs 
of these hopping terms coming from the eight-site clus- 
ter calculation are the same as in Eqs. (|) and (|). To 
compensate for neglecting three- and four-site terms, the 
hopping amplitudes in (0) are re-scaled by a single con- 
stant a and the amplitudes in (|^) by another constant f3. 
A reasonable fit of the energies of the effective model on 
the 4x2 cluster to the exact energies of the t-J model on 
the 4x4 cluster in the whole range studied, < J < 2.5, 
is achieved with a — 1 and /3 = 0.5. In this range of J, 
the relative difference between these two energies is less 
than 0.01. The exchange part of the effective model is 
given by Eq.(|). 

A first insight on the properties of the effective model 
can be gained by studying the 4x4 cluster. Due to 
the large dimension of the Hilbert space (« 7.76 10^ 
for Np ^ 1 and « 1.50 10^ for A^p = 2 taking into 
account translational invariance), conventional exact di- 
agonalization techniques cannot be applied, except by 
resorting to massive computers. Alternatively, a diag- 
onalization in a systematically expanded Hilbert space 
(SEHS)Ea is used. With a number of states « 5 10^ 
variational energies within 1% of the exact energies, es- 
timated by extrapolating to the full dimension of the 
Hilbert space, are obtained. In order to study larger 
clusters, a quantum Monte Carlo(QMC) technique wi^ 
the conventional worldline checkerboard decompositioncil 
is used. Although there are no fermions involved, there 
is a "minus sign problem" which makes impossible the 
study at low temperatures. There are several terms in 
the effective Hamiltonian which lead to "minus sign" con- 
figurations in the 2+1-dimensional space. An important 
reduction of this problem is achieved by not generating 
those configurations with interacting cubes which do not 
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conserve the parity of the number of triplets on its top 
and bottom plaquettes. Of course, the QMC algorithm 
is no longer exact but nevertheless it provides a reason- 
able approximation to its exact behavior. By eliminating 
certain transitions, it might be possible that certain re- 
gions of the phase space are disconnected. To cope with 
this problem, at each temperature, results were averaged 
over at least four independent runs starting from differ- 
ent initial states. Variations in the values of the energy 
from different runs were somewhat larger than the sta- 
tistical error of each run but nevertheless smaller than 
1% in all cases. In addition, since the the present study 
concerns zero temperature properties, the simulations are 
restricted to the subspace of zero total magnetization and 
only local moves are included in the algorithm. Overall 
simulation error bars are approximately twice the size of 
the symbols used. 
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FIG. 8. (a) Variation of energy per site with temperature 
in the 5*^ = subspace in the 4x4 cluster with one and two 
pairs, J jt = 0.8. (b) Variation of triplet density and (tTjTt) 
magnetic structure factor with temperature in the subspace 
5^ = in the 4x4 cluster with two pairs, J/t = 0.8. In both 
panels, solid symbols at T = indicate the corresponding 
values obtained by SEHS diagonalization. 

Evolution of various quantities with temperature in the 
S"^ = subspace on the 4x4 cluster, J = 0.8i, with one 
and two pairs are shown in Fig. ^. The corresponding 
zero temperature results obtained by diagonalization in 
an expanded Hilbert space are also included for com- 
parison. As expected the energies obtained by the ap- 
proximated QMC technique are higher than the obtained 
by diagonalization (which are virtually exact), while the 
triplet densities are smaller than the exact ones. 

The same correlations previously studied on ladders 
are shown in Fig. ^(a) for the 4x4 cluster, with one and 
two pairs, J — 0.8t. These ground state correlations, 
obtained by SEHS at k = (0,0), have been averaged 
over the two-directions x and y. The most important 



result is again that the triplets try to locate as far as 
possible from the pairs. Similar results were obtained for 
the other correlations shown in Fig. |^, and for all values 
of J jt examined. 

Results on the 8x8 cluster obtained by QMC tech- 
nique at T = 0.09t for one and two pairs are shown in 
Fig. ||(b). The behavior of triplet-triplet and triplet-pair 
correlations and pair-pair corrrelations are qualitatively 
similar to those found in the smaller cluster at zero tem- 
perature. These results are very suggestive that triplets 
are expelled from the vicinity of a pair. Notice also that, 
as for ladders (Fig. ||), singlet-singlet correlations decay 
very rapidly to their bulk value. 
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FIG. 9. Triplet-triplet (full symbols) and triplet-pair (open 
symbols) correlations as a function of distance, A'j, — 1,2, 
J = 0.8t, (a) on the 4 x 4 cluster (SEHS) and (b) on the 8 x 8 
one (QMC, T — 0.09t). Fair-pair (stars) and singlet-singlet 
(crosses, (b) only) correlations for A'j, = 2 are also shown. 

A more systematic study as a function of J/t leads to 
the results depicted in Fig. ^(a). In all cases, the prob- 
ability of finding a triplet near a pair (open symbols) is 
smaller than the average probability of finding a triplet 
on a given site (full symbols) . Consistently with the idea 
of pairs expelling triplet excitations, the probability of 
finding a triplet at the largest distance of a pair (only 
shown for the one pair case) is larger than the average 
triplet probability. Singlet-singlet correlations vs. dis- 
tance 

Since the system contains charged pair bosons, it is im- 
portant to calculate the superfluid density as a measure 
of superconducting properties in this model. To calculate 
the superfluid density it The procedure to calculate the 
superfluid densitynJ starts by computing the correlation: 



C{t) = (Ap,(T)Ap,(0)) 



(9) 



where Ap^ir) = Px{t +1) -p^ir) andpxiT) = J2^iii'^)i 
where x(i,r) is the ^-coordinate of pair i at imaginary 
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time T and the sum extends over all the pairs in the 
system. The superfluid density follows from: 



lirn I dTC{T)exp{—iujT) 



(10) 



It is well-known that this quantity vanishes if the simula- 
tion is carried in a zero- winding number subspace which 
is the case in the present study. The way out of this 
problem stems from the fact that the winding number 
may be nonzero in half of the Trotter, imaginary time, 
direction, being also nonzero but with an opposite sign in 
the other half. The superfluid density then would come 
out by taking the Fourier transform in (10) between zero 
and 13/2. This procedure should be exact in the limit 
of infinite Trotter number which is the limit in which 
on the other hand the whole worldline algorithm is valid. 
This procedure has been thoroughly checked in the hard- 
core bosoBiHipdel where exact and numerical results are 
available.Eay 
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FIG. 10. Results obtained by QMC on the 8 x 8 cluster, 
T = 0.09t at various fillings and as a function of J/t. (a) 
Triplet density (full symbols) and triplet-pair correlations at 
NN sites (open symbols). The triplet-pair correlation at the 
maximum distance for Np = 1 are shown with stars(b) AF 
structure factor (triangles: Np = 4, diamonds: Np = 8) and 
superfluid density (plus: Np — 4, stars: Np — 8). 

The superfluid density, ps, and the magnetic structure 
factor at {tt,tt), xaf, on the 8x8 cluster with four and 
eight pairs, are shown in Fig. [lo|(b) as a function of J/t. 
Although the ultimate dominance of SC or AF, or even- 
tually the coexistence of both SC and AF orders can be 
determined only by a finite size extrapolation which is 
out of the scope of the present study, two qualitative fea- 
tures are apparent. First, as a pair doping is increased 
there is an enhancement of SC and a suppression of AF. 
Second, as J/t is increased xaf also increases while p 
decreases. Similar results were obtained with an even 



simplified-,model which mimics the spin triplets by spin 
doubletsLJ and can be obtained by projecting out the 
S"^ = component of the triplets, tg, again by using 
(^). The mutually exclusive behavior of AF and super- 
conductivity, which has been shown both in the present 
work and previously in Rcf. can now be understood 
in terms of this real space phase separation between SC- 
singlet and AF regions. The growth of one of each phases 
reduces the space available for the other phase. 



IV. CONCLUSIONS 

In the first place, the present work is concerned with 
the mapping of the t-J model into an effective model of 
pairs and triplets moving on a "sea" of singlets. The 
strategy adopted was first to map exactly the magnetic 
interactions in the undoped system and then to deter- 
mine the effective hopping interactions in the two-hole 
sector of the original Hamiltonian. This second proce- 
dure is an approximate one because it implies projecting 
out single-occupied dimers and truncating the range of 
the interactions nearest neighbors. In the case of ladders 
this procedure is reasonably clean. The choice of dimers 
corresponds to the strong coupling limit which extends 
virtually to the anisotropic case. An effective model in 
one dimension is obtained in which the hopping couplings 
between pairs and singlets are larger than between pairs 
and triplets. 

In the square lattice the procedure to obtain the ef- 
fective model is more complicated. In the first place, 
the simple dimer covering adopted breaks rotational in- 
variance of the lattice. Still, the mapping of the mag- 
netic interactions at half- filling is exact. In the second 
place, even by restricting the range of hoppings to nearest 
neighbors, three- and four-site interactions appear at the 
effective level. To eliminate these interactions another 
fit to the energies of a small t-J cluster was performed 
to rescale the effective hopping constants. The resulting 
effective hopping couplings show the same behavior as 
those obtained for the ladder case. 

The resulting effective models were studied by numer- 
ical techniques, exact diagonalization in the case of the 
model obtained for the ladder lattice and approximate di- 
agonalization and Quantum Monte Carlo for the square 
lattice case. Again the results obtained in the effective 
model for ladders are more clear than the ones obtained 
for the square lattice. 

From the study of several correlations the following 
picture emerges. Pairs are surrounded predominantly by 
singlets, and triplet excitations are located as far apart 
as possible. This is a kind of phase separation between a 
pair-rich "RVB" phase and a pair-poor triplet-rich phase 
which would correspond to a phase with at least short- 
range interactions (since_triplet excitations restore AF 
order from a RVB statellj) and would quite likely be 
insulating. Pair-pair correlation functions indicate that 
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pairs try to be situated at the maximum distance (this 
behavior is more definite in ladders than in 2D). This 
behavior might correspond to pairs moving in a singlet 
phase macroscopically separated from the triplet phase. 
Alternatively, the phase separation might be microscopic: 
pair-singlet islands moving in a triplet background. The 
behavior of triplet-pair correlations in ID could favor this 
second scenario. In any case, these PS scenarios provide 
an explanation for the AF-SC mutual exclusion shown in 
the previous Section. 

Further studies are necessary to distinguish between 
both types of phase separation. In the macroscopic PS 
scenario it is expected a superfiuid density comparable 
with that of the hard-core boson model while on the mi- 
croscopic PS scenario it would be considerable reduced. 
On theoretical grounds, the emerging PS picture could 
be realted to the AF-SC coexistence phase predicted by 
SO (5) theories. The exclusion of triplets and pairs is 
highly nontrivial taking into account that previous stud- 
ies on the t-J modelcj have suggested a bound state be- 
tween a (Ix2_y2 pair and a triplet. However, the internal 
structure of a pair, essential in the analysis of Ref. |2^, 
is lost in the present study.E3 Finally, it is also tempting 
to relate this PS state to recent observations of inhomo- 
geneities in Bi2Sr2CaCu208-i-5 (Ref. [l^. It might be pos- 
sible that out-of-plane negative Coulomb centers could 
attract and pin pair-singlet islands in their surroundings. 
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